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1. Introduction 


Meson masses are not only fundamental quantities of hadrons but also a key to know properties 
of quantum chromodynamics (QCD) vacuum. At finite temperature (T), we can define fwo kinds 
of meson masses, pole and screening mass. Meson pole masses are one of fhe possible observables 
in fhe heavy ion collisions. Screening masses of lighf mesons are essenfial for fhe range of fhe 
nuclear force. Accordingly, if is necessary fo consfrucf fhe effecfive model for calculafing pole and 
screening mass simulfaneously. 

In lattice QCD(LQCD), meson pole (screening) masses are calculafed from fhe exponenfial 
decay of femporal (spatial) mesonic correlafion funclions. LQCD simulations are more difficull for 
pole masses fhan for screening masses, since fhe lattice size is smaller in fhe fime direcfion fhan in 
fhe spatial direcfion. This sifuafion becomes more serious as T increases. For ibis reason, meson 
screening masses were calculafed in mosf of fhe LQCD simulations. Recenfly, a slale-of-lhe-arf 
calculation was done for meson screening masses in a wide range of T < 800 MeV |[T]] 

Consfrucfing fhe effecfive model is an approach complemenfary fo fhe firsf-principle LQCD 
simulation. In confrasf fo LQCD simulafions, meson pole masses are exfensively invesfigafed af 
finite T by fhe Nambu-Jona-Lasinio (NIL) model |Q, fhe Polyakov-loop extended Nambu-Jona- 
Lasinio (PNJL) model [Q]. However, only a few frials were made so far for fhe evaluation of meson 
screening masses here ^ means a species of mesons. The model calculations have 

essenfially fwo problems. One problem is fhaf fhe NJL-fype models are nonrenormalizable and 
hence fhe regularization is needed in fhe model calculations. The regularization commonly used is 
fhe fhree-dimensional momenfum cutoff. The momenfum culoff breaks Lorenfz and franslafional 
invariance, fhereby fhe spatial correlation funclion T]^^(r) has an unphysical oscillafion [^. This 
makes fhe determination of quite difficulf, since M^,scr is defined from fhe exponenfial decay 
of (r) af large disfance (r): 


= - lim 

r—>£>o 


‘ailnT7jg(r) 

dr 


( 1 . 1 ) 


Anofher problem is fhe feasibilify of numerical calculafions. In fhe model approach, (r) is 
firsf obfained in fhe momenfum (q) represenfafion lu fhe Fourier fransformafion fo fhe 

coordinafe represenfafion (r = |x|). 


nuir) = / myo.f V'” . ( 1 . 2 ) 

The infegrand is slowly damping and highly oscillafing particularly af large r where is de¬ 

fined. This requires heavy numerical calculafions. If was fhen proposed fhaf fhe contour integral 
was made in fhe complex-^ plane [^. However, fhe confour integral is sfill hard fo do because of 
fhe presence of fhe femperafure cufs in fhe vicinify of fhe real axis [||] ; see fhe leff panel of Fig. |^, 
where nofe fhaf e is an infinifesimal quantify. 

In fhis falk, we propose a new formalism for calculafing screening mass and discuss fhe pos- 
sibilify of fhe predicfion for pole mass from screening mass by using effecfive model. This falk is 
based on fhe paper |Q]. 
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Fig. 1 : Singularities of in the complex-^ plane based on the previous formulation [|^] (left) and 

the present formulation (right). Cuts are denoted by the wavy lines and poles by the points. 

2. Formalism 

The Lagrangian density of the two-flavor EPNJL model ||^ is defined as 

^ = q{iYvD'' - mo)q + Gs{<P)[iqq)^ + {qiYs^qf] - ^i<P[A],0[A],T) (2.1) 

with the quark field q, the current quark mass mo and the isospin matrix f. The coupling constant 
Gs(<P) of the four-quark interaction depends on the Polyakov loop 0 as 

Gsi0) = G, [l-ai^'^-a2(^'^ + ^^)] , (2.2) 

where = d^' + iA^ with = 5og{A^)aK/^ = —5Qig{A4.)aXa/2 for the gauge field A^, the 
Gell-Mann matrix Xa and the gauge coupling g. When ai = a 2 = 0, the EPNJL model is reduced 
to the PNJL model [Q]. 

In the EPNJL model, only the time component of A^ is treated as a homogeneous and static 
background field, which is governed by the Polyakov-loop potential The Polyakov loop 0 and 
its conjugate 0 are then obtained in the Polyakov gauge by 

^=itr,(L), a = itrc(L*) (2.3) 

with L = exp[;A 4 /r] = exp[/diag(A{\A|^,A 4 ^)/r] for the classical variables A'l satisfying that 
A^ +A^ -h A 4 ^ = 0. Eor zero chemical potential, 0 equals to 0. Hence it is possible to set Ap = 0 
and determine the others as A^ = —A 4 ' = cos^'(2^^)r. We use the logarithm-type Polyakov- 
loop potential 'fX of Ref. [^, but refit the parameter Tq to reproduce the chiral phase transition 
temperature Tc because the original value of To is set to 270 MeV which is the deconfinement 
transition temperature in the pure gauge limit. 

Making the mean field approximation to ( |2.1[ ) and the path integral over the quark field, one 
can get the thermodynamic potential (per unit volume) as 

[1 -h 3(^> -h 

[1 -h 3(^ -h 

(2.4) 


Q =Um + '^-2Nf 


/ 


d^p 

( 271)3 


3E„ 


+ jln 


+ jln 
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with j8 = 1 /T, M = mo — 2Gs(^*)c7, Ep = s/iP-+ M^, and Um = Gs{<P)(y^. Here, a means chiral 
condensate {qq). Nf is the number of flavors. We determine the mean field variables {X = a, 0, <P) 
from the stationary conditions for 



(2.5) 


Since the momentum integral of diverges, we use the Pauli-Villars (PV) regularization [^, 
In the scheme, the integral I{M,q) is regularized as 




a=0 


( 2 . 6 ) 


where Mq = M and Ma (a > 1) are masses of auxiliary particles. The parameters Ma and Ca should 
satisfy the condition La=o^ce = La=o^ce^a = 0- We then assume (Co,Ci,C 2 ) = (1,1,—2) and 
(MpMl) = (M^ + 2A^,M^ +A^). We keep the parameter A finite even after the subtraction (^!^), 
since the present model is nonrenormalizable. The parameters taken are mo = 6.3 MeV, Gg = 5.0 
GeV^^ and A = 0.768 GeV. This parameter set reproduces the pion decay constant fji = 93.3 MeV 
and the pion mass Mj^ = 138 MeV at vacuum. 

We derive the equations for pion and sigma-meson masses, following Ref [Q]. We consider 
currents with the same quantum number as pion (P) and sigma-meson (S), 


7p“(v) = q{x)iy 5 ffq{x) , Js{x) = q{x)q{x) - {q{x)q{x)). 


(2.7) 


The Fourier transform of the mesonic correlation 


function T]^^(v) = (0|r \ |0) i 


IS 


= i J d\ P‘?-"(0|T (7^(v)7^^(0)) |0), (2.8) 

where ^ =P‘‘ for pion and S for sigma meson and T stands for the time-ordered product. Using the 
random-phase (ring) approximation, one can obtain as follows, 

^ i-2Gg(^)n^^’ 

where the one-loop polarization function is explicitly obtained by 

Hss = 2iNf[h +l2- {q^ - 4M^)h] , Hpp = 2iNf[Ii +I 2 - q^h], (2.10) 


with 


h = 
h = 


tTr 


1 


d'^p 

(271)4 “"Lp'2_ m2 J 

d^p 


rtTr 


> h = 

1 


d^^p 

(271)4 L (p' -I- ^)2 — m2 . 


trn 


(271)4 n{(p'+^)2-M2}(p'2_M2)J’ 


( 2 . 11 ) 

( 2 . 12 ) 


Here, q^ = q^ — q^ and p' = {po -h iA^,p). ttc means the trace of color matrix. For finite T, the 
corresponding equations are obtained by the replacement 


Po —^ i(Oi — i(2l l)7rP , 



00 

—> iT ^ 

I —— 00 


d^p 

(W 


(2.13) 
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The meson pole mass p^jg is a pole of Taking the rest frame q = (^O)O) for 

convenience, one can get the equation for p^jg as 


[l-2Gs(^')n^^(^^,0)]| 


= 0 . 


(2.14) 


I go ,pole 

The method of calculating meson pole masses is well established in the PNJL model |0]. 

The meson screening mass gg^ defined with ( |i~T| ) is obtained by making the Fourier trans¬ 
form of shown in ( |1.2[ ). In the previous formalism ||3|], however, the procedure re¬ 

quires heavy numerical calculations in the part, as shown below, where means a function 
after the PV regularization. Taking the I summation before the p integral in ( 2.13 ), one can describe 
as the sum of the vacuum and temperature parts, and , defined by 




InM^ +/vi 


2Ma 


(0 


£ C« rd\p\ ftU\p\,q) [F+{Ep)+F-{Ep)] , 

a=0 


/vac(-r) = V 1 +X^ln 


Vl +X^ + 1 


, /tem(|p|,^) = 

Ep q 


Vl +x^ — 1 

where E± are the Fermi distribution functions. E± are defined as 


(q-mf + e^ 
(q + HpIV + e^ 


(2.15) 

(2.16) 
(2.17) 






1 


(2.18) 


Nc g + 1 

integral well defined at ^ = ±2|p|, but this requires 


In ( 2.17 ), the term is added to make the 
the limit of £ —)• 0. 

As shown in the left panel of Fig. |T|, ffacV-Ma/ q) and /tem(|p|,^) have the vacuum and tem¬ 
perature cuts in the complex q plane, respectively. In (pT^, the cuts contribute to the q integral in 
addition to the pole at ^ defined by 


[l-2Gs(^>)n^^(0,^2 




= 0 . 


(2.19) 


It is not easy to evaluate the temperature-cut contribution, since in ( |1.2[ ) the integrand is slowly 
damping and highly oscillating with q near the real axis in the complex q plane. Furthermore we 
have to take the limit of £ —>• 0 finally. In order to avoid this problem, we integrate about p in (2.13) 
before taking Matsubara summation Consequently, we can rewrite as an infinite series of 
analytic function. 




iT 

4nq 


Nc oo 2 

L L £Casin^^ 


!=1 l=~<x>a=0 



( 2 . 20 ) 


where 

Mu,a{T) = ^MI + {{21 + 1 )kT+A^F (2.21) 

We have numerically checked that the convergence of I summation is quite fast in ( [2.201 ). Each 
term of has only two cuts starting from EliMipa on the imaginary axis in the complex q 
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plane. The cuts are shown in the right panel of Fig. |I|. The lowest branch point is ^ /=o,a=0- 

Hence 2M,=i /=o,a=o is regarded as “threshold mass” in the sense that the meson screening-mass 
spectrum becomes continuous above the point. 

If scr < 2M^■=I,/=o,a=o^ the pole at ^ is well isolated from the cut. Hence one 

can take the contour (A—;-B—>-C—)-D—)-A) shown in the right panel of Fig. [T]. The q integral of 
{0,q^)e'^'' on the real axis in ( [T^ ) is then obtained from the residue at the pole and the line 
integral from point C to point D. The former behaves as exp[—M^^^^^r]/r at large r and the latter as 
exp[—2M,=i/=o,a=o^]/'"- The behavior of (r) at large r is thus determined by the pole. One can 
then determine the screening mass from the location of the pole in the complex-^ plane without 
making the q integral. In the high-T limit, the condition tends to < 2nT. 

3. Numerical Results 

The pion screening mass obtained by state-of-the-art 2 +1 flavor LQCD simulations [|^ 
is now analyzed by the present two-flavor EPNJL model simply, since pion is composed of u and 
d quarks. This is a quantitative analysis, because the finite lattice-spacing effect is not negligible 
in the simulations. The chiral transition temperature is evaluated as = 196 MeV in the sim¬ 
ulations [p, although it becomes = 154±9 MeV in finer 2-1-1-flavor LQCD simulations 
close to the continuum limit. Therefore, we rescale the LQCD results of Ref. [|I[] with multiplying 
them by the factor 154/196 to reproduce Tc = 154 ±9 MeV. The model parameters, mo and Tq, 
are refitted to reproduce the rescaled 2-1-1 flavor LQCD data, i.e., Mji = 175 MeV at vacuum and 
Tc = 154 ±9 MeV; the resulting values are mo = 10.3 MeV and To = 156 MeV. The variation of 
mo from the original value 6.3 to 10.3 MeV little changes a and 

As shown in Lig. the M^^.scr calculated with the EPNJL model (solid line) well reproduces 
the LQCD result (open circles), when a\ = a 2 = 0.31. In the PNJL model with ai = a 2 = 0, the 
model result (dotted line) largely underestimates the LQCD result, indicating that the entanglement 
is important. The dashed line denotes the sigma-meson screening mass scr obtained by the 
EPNJL model with a\ = a 2 = 0.31. The solid and dashed lines are lower than the threshold mass 
2M,=i /=o,a=o (dot-dashed line). This guarantees that the M^^ scr and determined from the 

location of the single pole in the complex-^ plane agree with those from the exponential decay of 
at large r. The chiral restoration takes place at T = = 154 MeV, since M^r.scr = 4/(j,scr 

there. After the restoration, the screening masses rapidly approach the threshold mass and finally 
27iT . The threshold mass is thus an important concept to understand T dependence of screening 
masses. 

Linally, we predict the T dependence of pole mass p^ig for pion and sigma-meson with 
EPNJL model (Lig. At low temperature {T < Tc), the T dependence of p^ig and are 

almost same in the pion and sigma-meson because Lorentz symmetry is preserved approximately. 
Around Tc, pion and sigma-meson masses agree with each other and chiral symmetry restoration 
takes place at the same temperature Tc for pole and screening mass. These indicate that at low tem¬ 
perature {T < Tc) we can predict the T dependence of p^jg from that of simply. Above Tc, 
however, the difference between p^jg and ^g^ gets larger as temperature increases. Therefore, 
above Tc, it is necessary that we should use the effective model to predict the pole mass from the 
lattice QCD results of screening mass. 
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Fig. 2: T dependence of pion and sigma-meson Fig. 3: T dependence of screening and pole mass in 
screening masses, and M^^scr- pion and sigma-meson. 

4. Summary 

We have proposed a practical way of calculating meson screening masses in the NJL- 
type models. This method based on the PV regularization solves the well-known difficulty that the 
evaluation of is iioi £^sy in the NJL-type effective models. In the previous formalism |Q], 

the vacuum and temperature cuts appear in the complex-^ plane. The contributions to the mesonic 
correlation function are partially canceled in the present formalism. The branch point of the re¬ 
sultant cut can be regarded as the threshold mass. The pion and sigma-meson screening masses 
rapidly approach the threshold mass 2M,=i,/=o,a=o(^) chiral restoration. We propose the 

prediction for pole mass from screening mass by using EPNJL model. 
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